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We show how to implement quantum computation on a system with an intrinsic Hamiltonian by 
controlhng a Hmited subset of spins. Our primary result is an efficient control sequence on a nearest- 
neighbor XY spin chain through control of a single site and its interaction with its neighbor. Control 
of an array of sites yields sufficient parallelism for the implementation of fault-tolerant circuits. The 
framework exposes contradictions between the control theoretic concept of controllability with the 
ability of a system to perform quantum computation. 



o 



Oh 
I 



> 
o 

o 

in 
o 

On 
O 

> 

X 



Introduction: What does it take to implement a quan- 
tum computation in a given physical system? This would 
seem to be a fundamental question, for which a sufficient 
set of conditions is well known [Tj ; implementation of sin- 
gle qubit rotations on any spin, and a nearest-neighbor 
two-qubit gate. However, since that degree of control 
seems to be extremely demanding, it is vital to under- 
stand how little control is required. In fact, there are 
some systems whose internal dynamics are sufficient to 
implement computations [2J[3]. However, these have to 
be carefully designed, and still require the ability to pre- 
pare the initial (product) state. On the other hand, rein- 
troducing control over a single spin in principle gives suf- 
ficient control for almost all Hamiltonians 4 . This archi- 
tecture is an attractive proposition for some experimental 
implementations. For instance, while it has been shown 
that full control over a pair of superconducting qubits 
can be achieved, the physical layout of such devices is 
highly asymmetric [S], and is not easily scaled up. Thus, 
we might be able to consider manufacturing a uniform 
system with some fixed interaction, and can concentrate 
all our design efforts in generating controllability at the 
end of a chain, which is allowed to be non-symmetric. 



The proofs of controllability of these interface schemes 
[4] make no claims regarding efficiency. Some examples of 
Hamiltonians have been specifically constructed to allow 
efficient control sequences [21 Ej . While much less com- 
plicated than those of [H [3] which function without any 
control, they are still unrealistic. In this paper, we de- 
velop efficient analytic control sequences for a much more 
natural class of Hamiltonians; spin chains. The main in- 
gredients are an encoding of information in the diagonal 
basis of the Hamiltonian, and the use of Rabi oscillations 
to induce transitions between these states. 



Generic Controllability: Consider an A'^-qubit Hamil- 
tonian H, with control field hi. Each arbitrarily entan- 
gled eigenvector \Xx) can be identified with a logical basis 
state \xl), X € {0, 1}^. Generically, the eigenvalues |Aj;| 
and differences [Ax — Aj,| are unique, and {X^l hi \Xy) ^ 0. 



Under these assumptions, the field 
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applies the logical X rotation on qubit n (up to some 
phases, which we consider later). x®nis used to denote 
the flipping of bit n in the string x. Naturally, h\ only 
makes sense if {\x®n\ I'^k) 7^ 0- Due to the assumed 
uniqueness of gaps, each term in the sum is on resonance 
with a single transition so that, by applying the rotating 
wave approximation (RWA, which requires that the de- 
tunings of different energy gaps is much greater than B), 
the effective Hamiltonian in the interaction picture is 



Hrff — B 



>(A. 



which evidently provides the logical X rotation that we 
desire, by any angle Btx- The effect of returning to the 
Schrodinger picture is that this rotation is followed up by 
^^g-»A,tx |;^^^ ^;^^|^ A cNOT gate (up to an identical 
phase condition) is implemented in a similar fashion. 
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with control qubit m and target n. In order to have full 
controllability, we just need to demonstrate how to imple- 
ment arbitrary Z rotations on any spin, n. This can also 
be used to cancel the phases that accrue due to the inter- 
action picture. The first step is to negate the effect of the 
phases when implementing an identity operation, using 
the standard NMR technique of refocusing - by applying 
the cyclic permutation |A^_|.i mod 2«) {^x\ 2^ times, 
waiting the same time tz between each application, then 
all eigenvectors accumulate the same phase, tzTr{H). In- 
dependently varying the waiting times in different inter- 
mediate states allows different phases to be applied to 
different eigenvectors, which is precisely what we need, 
thereby proving controllability of a generic Hamiltonian. 

This technique is, in the majority of cases, wildly in- 
efficient, for several reasons. Primarily, since there is an 
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exponential number of eigenvectors, B must be exponen- 
tially small if the control field is to be bounded, so gates 
take exponentially long. Equally, to cycle through all 
the eigenvectors for the phase gate is an exponential pro- 
cess. While these techniques are not necessarily unique, 
the on-resonant control would seem to be an essential 
component of any such scheme. How can any scheme be 
efficient? Introducing degeneracies into the system re- 
duces the number of terms that we sum over. However, 
care is required since, if we have that A^^ — Xx(Bn is inde- 
pendent of X (such that is only a single term), then 
there is too much degeneracy for hcNOT, and the existing 
proof of controllability breaks down. Bizarrely, to get ef- 
ficient computation, we have to make it harder to prove 
controllability! Worthy of emphasis is that controllabil- 
ity typically applies to the control of the entire Hilbert 
space, whereas efficient quantum computation only re- 
quires control over a subsystem. 

Computation on spin chains: While it may be inter- 
esting to understand generic systems, the Hamiltonians 
that are accessible in the laboratory are far from generic, 
so the preceding arguments need not apply. We shall now 
show how to use the basic ideas introduced to efficiently 
compute on a spin chain of the form 

N N 
71—1 n—1 

(1) 

This Hamiltonian is exactly solvable [71, the first step 
being to perform the Jordan- Wigner transformation aj^ = 
'^n Y[m=i ^m- Thcsc Can then be transformed into a set 
of non-interacting fcrmions, 

N 
n=l 

via a Bogoliubov transformation and diagonalization of 
a.n N X N tridiagonal matrix. We shall assume that the 
coupling strengths J„ are known, although they can be 
identified experimentally [8] after preparing the system 
in some initial state [9j. For pedagogical reasons, we 
introduce three control fields, which only act on the first 
two spins [T5], 

hi - Xi 

h2 = i((l+7)^X + (l-7)yy)i,2 

/l3 = ZiZ2- 

These can be transformed into the {a„} or {6„} basis. We 
can assume Ji = Bi — 0, which means that bi = oi and 
Ai = 0. In cases such as 7 = 0, it is already guaranteed 
that the {A„} are unique and that Q!„ := (0| 6i/i2fej^ |0) 7^ 
[To], where |0) denotes the vacuum (i.e. ground) state 
of the system. The uniqueness of the {A„} is sufficient 
to give the condition of uniqueness of the {|A„|} since we 
could tune the field Bi which, working in an offset system 



where we keep Ai =0, rescales all other eigenvalues by 
Bi, sufficient to move them off any degeneracies due to 
the existence of ±A„ eigenvalue pairs. It is also sufficient 
to ensure that none of the eigenvalues are exponentially 
small. Henceforth, we assume these conditions hold. 

Instead of proving universal computation on the full 
Hilbert space, we shall just consider a subspace where 
the logical qubits are described by pairs of fermions. The 
initial state is of the form 

L(7V-1)/2J 

io,^«L(--i)/2j ^ n 4jo), 

and the raising operator for the n*^ logical qubit is 
CT+ — ^'2n+i^2n- The primary reason for this choice is that 
if we were to encode in single fermion states, then when 
moving states around the lattice, they generate exchange 
phases, which correspond to controUed-phase gates. En- 
coding in a |01)^,|fO)^ subspace negates these effects 
[llj . Note that bi is not used to encode a qubit, and is 
instead kept free, as workspace. 

All protocols in the computation require the field 

B.n{t) = Bcos(A„i)/i2, 

which implements the effective Hamiltonian 

^anB{b\bn + bibi). 

The corresponding unitary evolution is just a swap be- 
tween the two modes 61 and 6„, except that when we 
reassert the normal ordering of the fermionic modes, a 
phase factor of HrniiC^^L^m ~ 1) arises if the swap oc- 
curred. The sequence of 111^=2(2^™^™ ~ ^) precisely 
the c-phase gates mentioned previously, whose effects are 
negated by the encoding - that term calculates the parity 
of the number of fermions in modes 2 to n — 1 if there's 
only 1 fermion in modes 1 or n, and this number is fixed 
due to our encoding. Thus, up to a diagonal gate, Bn(t) 
can be used to implement a swap of a fermion in mode 
n onto spin 1. When implementing this swap, one of the 
two states will always be empty, so the diagonal gate is 
only a local phase gate, which we will later see how to 
correct (either we swap a fermion onto the empty state 
on site 1, or we undo that swap). Once we have imple- 
mented B2n{'^ / {Ba2n)) to swap fermion 2n to the first 
site, we can implement i?2n+i(26'/(i?a2n+i)) before ap- 
plying B2n{T^ / {Ba2n)) ■ This returns the fermions to their 
original positions but the logical qubit n has undergone 
an X-rotation of angle 9, up to the phase gates due to the 
transformation between the interaction and Schrodinger 
pictures. This protocol also allows the preparation of any 
eigenstate of Hf and measurement of any logical qubit; 
swapping n to 1, measuring and swapping back projects 
the system into a Fock state of 6„, and hi allows the bit 
to be flipped after measurement. Fig.[l]demonstrates the 
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FIG. 1: The effectiveness of the swap gate between the bi 
mode and the 6„ with minimum eigenvalue for A'^ — 101 and 
the couphng scheme of Eq. ^ , offset by using Bi = \/2. The 
inset shows how decreasing B increases the accuracy. The 
fidelity is calculated as the overlap between the evolution of 
an initial state |0) |Aioo) and the target state |1) |0)®^°". AU 
quantities are dimensionless by taking h — 1. 

simple swapping protocol for a chain of 101 spins in the 
single fermion subspace. 

A refocusing technique can now be used to perform ar- 
bitrary Z rotations. If an X rotation is performed on each 
logical qubit every tz, then they each acquire a global 
phase of the form {X2n + •^2n+i)'7it2 a-t times 2mtz- By 
performing the gate X„ at times t'^ (instead of tz) and 
2tz, we get a phase rotation of 2(A2n+i — \2n){tz ~ i'z)- 

To entangle logical qubits m and n, we apply 
S2m(7r/(_Ba2m)), Swapping the fermionic mode 2m onto 
spin 1, followed by a field B'J^l) = 2i?'cos((A2n — 
A2n+i)0^3 fo'' time 9/{2B')^ which gives an effective in- 
teraction between modes 2n and 2?! + 1 (in the interac- 
tion picture), dependent on the presence or absence of a 
fermion on the first spin, 



ifeff = B' a2nOl2n+l{'^b\bi - l){h\J)2n^ 



4n+1^2ri) 



By applying B2m{T^ / {Ba2m)) , the sequence is completed. 
The ultimate result is a c-X rotation of angle 9, targeting 
qubit n, up to local rotations. 

This proves the possibility of implementing computa- 
tional gates on a sufiiciently large subspace. However, it 
is not sufficient for efficiency since the timing condition is 
based on the requirement that B and B' are sufficiently 
small. This gives two conditions to satisfy, B < \ and 
B max q;„ <C min | A„ — Am | ■ The first of these arises from 
the desire to only use finite field strengths. If any eigen- 
values, or their gaps, are exponentially small, or overlaps 
of eigenvectors on the second spin are exponentially small 
(any of which can happen, albeit rarely), then the gate 
time must be exponentially long. This loss of practi- 
cal controllability as a theoretically controllable system 
closely approaches a symmetric uncontrollable system 
has recently been identified in T^. In the case of the 
uniformly coupled chain (7 = 0, J„ = 1,-B„ = 0), the 
detunings are of the order of so gate times are 



introduced in [TO], with 7 = Bn = 0, 

3n2((Ar-l)2- 



/2 



N{N -2){2n- l)(2n + l) 



(2) 



It has a spectrum with regular spacings of 2/(iV — 2) 
and a„ = X/y/N — 1, meaning that it can implement 
gates in a time 0(N), which is optimal if J ^ 0(1). 
Some care has to be taken with the two-qubit gate, since 
the gaps between eigenvalues are highly degenerate. The 
first step in overcoming this is to make a suitable associ- 
ation between the numbering of the fermionic modes and 
their eigenvalues, X2n = — 1 + 2(n — l)/(iV — 2), A2n+i = 
{2n —1)/{N — 2), which means that applying B'^{t) uses 
a frequency greater than half of the total energy range 
i.e. each mode can only couple to one other mode, and 
corresponds to applying the same X rotation on every 
logical qubit simultaneously, if the bi mode is occupied. 
To localixe this effect on a single target, we apply a logical 
Z gate at the start, and half way through the evolution, 
on the qubits where we don't what the cNOT applied. 
Since ZVxzVx = 1, the evolution is canceled, and if 
the 61 mode was not occupied, we only get ZZ = 1. 

The engineered coupling scheme of Eq. ^ is partic- 
ularly amenable to the final step of the analysis - an 
estimation of how the gate error scales with TV and B. 
Instead of considering _B„ (t) , we will replace it with 



Bn{t) 



B 



cos(A„i)((i + i)xx + (1 - 7)rr)i 



+ sin(A„<)((l + 7)Xr - (1 - ^)YX)i,2, 



reducing our reliance on the RWA, 
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B 



2VN^ 



N 
m=2 



iA„t 



blbi- 



blb„ 



neglecting, for convenience, the string of operators 
J^^^2^(26^6fc — 1) which we know to be irrelevant due 
to our choice of encoding. A suitable rotating basis can 
be chosen to entirely remove the time dependence. As a 
first step, we estimate the error in the rotation within the 
{5i,5„} subspace by adiabatically eliminating the other 
levels. This leads to an error of e ~ B^N \og^ N, using 



N 

E- 



1 



^'^ N-2 
< — A^logiV. 

m— 1 



2m 



We also need to estimate the leakage out of this subspace, 
which can be achieved by assuming the desired evolution 
of the subspace, in particular the amplitude of the b\ 
mode can be taken to be cos{Bt/^/N). Using this, the 
evolution of the other modes can be solved exactly, and 
their maximum amplitudes can be bounded. Summing 
0{N'^). Superior schemes can be designed, such as that all these reveals a maximum error of £ - B^N log^ N. 



4 



Therefore, by selecting B ^ ^ \/Nlog{N), the error is 
held constant, and the gate time scales as 0{N log N). In 
Fig. [ij the error e is evaluated numerically for fixed N, 
and indicates e ^ B^'^. Other coupling schemes are more 
strongly affected, but a choice of B ~ min(a„) min | Am — 
A„|/iV ensures a constant error with increasing N. 

One might ask how robust this scheme is to fluctua- 
tions in the control fields. Since we are using Rabi oscil- 
lations, there is a lot of built-in tolerance - the pulse 
sequence can be anything provided it has the correct 
Fourier component with the correct amplitude. Other 
Fourier components are irrelevant provided they are suf- 
ficiently far from the other energy gaps of our system. If 
the (integrated) amplitude of our Fourier component is 
slightly wrong, then that means the angle of the imple- 
mented X rotation is incorrect by the same fraction, but 
it is exactly the same possibility of error that all non- 
topological schemes suffer from. Similarly, if there is a 
small (compared to Ban) frequency discrepancy, this in- 
troduces a small Z-component to the X rotation. The 
Fourier decomposition of control sequences also indicates 
a link with [T3] where control of the single excitation 
subspace was demonstrated. Evidently, our fields i3„(t) 
give efficient controllability of this subspace for any spin 
preserving network, via Givens rotations (the exchange 
phases never manifest in the single excitation subspace). 
In [T3], the numerical techniques which suggested effi- 
ciency were based on a simple on/off switching of /i2, 
which can be directly related to our result by examining 
the Fourier modes of the square wave. 

Fault tolerance: This interface scheme has many ad- 
vantages such as not needing to perfectly engineer the 
system to within tight constraints. Instead, system to- 
mography can feed back into the control sequences. Also, 
at least theoretically, the majority of the system can be 
isolated from the environment, thereby decreasing deco- 
herence. Nevertheless, the possibility of error correction 
remains a concern. This introduces a significant problem 
to the interface scheme; as the system size increases, the 
errors accumulate more rapidly than they can be cor- 
rected. However, the architecture described here read- 
ily generalizes to structures with sufficient parallelism 
for fault-tolerance [Ij. Consider the system of Eq. ([iJ, 
but where we control some fixed set of spins {h}, by 
which we mean that we control the spins ki and /12 and 
/13 couplings between neighboring pairs {ki — l,ki) and 
{ki,ki + 1). By considering the scenario where all these 
couplings are switched off, the basis defines the computa- 
tional basis. On each site, if we only ever allow one of the 
couplings {ki — 1, ki) or {ki, ki + 1) to be active at a time, 
gates can be implemented in time 0{1) within a block, or 
between neighboring blocks. This is sufficient to design 
a fault-tolerant scheme [14], although care is required 
since errors that occur independently on each physical 
qubit correspond to correlated errors in the encoded ba- 
sis, constrained within a specific block. The constant 



sized blocks can be arranged into any geometry, allowing 
improvements in the fault-tolerant threshold. One would 
expect a threshold for per spin error rates of the order of 
Sc/K where Sc is any fault-tolerant threshold constrained 
by a locality condition, and K is the number of spins in 
any given block. 

Conclusions: Simple systems of non-interacting 
fermions, which can be converted to a wide variety of 
spin models, including XX and transverse Ising, can 
be efficiently controlled through the coupling of a single 
spin to its neighbor, enabling implementation of quan- 
tum computation. Without the additional coupling, the 
structure of the Hilbert space is entirely described by rep- 
resentations of SU(iV), which can be simulated in poly- 
nomial time on a logarithmic number of qubits, but in- 
troduction of a single controlling interaction breaks this 
symmetry and potentially permits a computation. The 
remarkable aspect is the ability to present analytic, ef- 
ficient, pulse sequences to achieve a computation. We 
have further discussed how the result generalizes to an 
array of controllers, which are sufficient to allow a fault- 
tolerant implementation; a feature absent from previous 
constructions [3J |2] • Our formalism motivates the expec- 
tation that most systems, while controllable, cannot be 
efficiently manipulated. This includes many interesting 
systems such as Heisenberg chains. 

In parallel to this work, Burgarth et al. have consid- 
ered the same problem |15) . In essence, our work proves 
when good solutions exist, at which point jl5j can be 
used to numerically find control sequences with smaller 
overheads (no proof for the existence of, or efficient con- 
vergence to, solutions is given in [rSj). 
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